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Fast 3D Synthetic Aperture Radar Imaging from 
Polarization-Diverse Measurements 
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Abstract 

An innovative 3-D radar imaging technique is developed for fast and efficient identification and characterization 
of radar backscattering components of complex objects, when the collected scattered field is made of polarization- 
diverse measurements. In this context, all the polarimetric information seems irretrievably mixed. A direct model, 
derived from a simple but original extension of the widespread “multiple scattering model” leads to a high 
dimensional linear inverse problem. It is solved by a fast dedicated imaging algorithm that performs to determine at 
a time three huge 3-D scatterer maps which correspond to HH, VV and HV polarizations at emission and reception. 
It is applied successfully to various mock-ups and data sets collected from an accurate and dedicated 3D spherical 
experimental layout that provides concentric polarization-diverse RCS measurements. 

Index Terms 

3D radar imaging, polarization, multiple scattering model, linear inverse problem, regularization. 


I. Introduction 

S YNTHETIC Aperture Radar (SAR) imaging is a widespread technique to make a map or image of 
the spatial distribution of reflectivity or backscattering of an object or scene from measurements of 
the scattered electric held. It is of great importance in many areas: land mapping, target recognition, 
environment monitoring, surveillance, nondestructive testing, etc. Radar imaging can also be used in an 
analysis purpose. For instance, to reduce the interference of wind turbine blades with air traffic control 
radars Q or to characterize tree signatures for canopy monitoring Q. It is then called ’’Radar Cross 
Section analysis” or “RCS analysis”. Regarding the data, it is usually collected from a remote platform, 
for instance a satellite or an airborne. It can also be collected nearby, e.g. inside an indoor facility, if the 
object is small enough, for analysis imaging. Most of the time, radar imaging is one-dimensional (1-D) 
or two-dimensional (2-D), leading respectively to 1-D backscattering profiles or 2-D backscattering maps. 
Three-dimensional (3-D) radar images are known to be far more complicated and challenging to process 
0 beyond the required extensive signal processing. Indeed, if a two-dimensional (2D) backscattering 
image can be formed by synthesizing an 1-D aperture with a wide-band radar, a three-dimensional (3D) 
backscattering image needs an entire 2D aperture. More precisely, in an indoor anechoic chamber, typical 
2D aperture geometries are planar, spherical and cylindrical Q. Their physical implementations require 
many illumination viewpoints around the object. It requires to combine many object and/or antenna 
rotations, a time-consuming task while the measurement conditions may derive. Once the backscatter 
field data is recorded, it is processed in order to get 3D images of the target radar backscattering spatial 
distribution. 

Many techniques have been developed for 3D radar image formation, both in SAR where the radar 
platform is moving while the target stays motionless and ISAR (Inverse SAR) where the target is moving 
while the radar platform radar is stationary. High resolution (HR) radar images are commonly obtained 
by processing coherently the backscattered fields as a function of the frequency and the angle (i.e. object 
attitude relatively to the radar). Among 3D HR radar imaging techniques emphasized in Q, there are the 
Polar Format Algorithm (PFA), also known as Range-Doppler, the Range Migration Algorithm (RMA), 
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its 2D version being also known as k-oo algorithm, and the Chirp Scaling Algorithm. Based on the polar 
nature of the frequency-domain backscatter data, 3D PFA is extensively used. Taking advantage of the 
processing that reduces in far-field condition to a Fourier synthesis problem, it is practically achieved 
via an interpolation, that reformats the data in the spatial frequency domain, and an inverse fast Fourier 
transform (FFT). Refer to Q or Q for implementation details on 3D PFA, including the polar reformatting 
mapping technique. Regarding RMA, it comes from seismic engineering and geophysics. Based on 1-D 
Stolt interpolation and again FFT through an approximation with the method of stationary phase (MSP), 
RMA is able to compensate completely a potential wavefront curvature. The last one, i.e. the Chirp Scaling 
Algorithm, is widely used in airborne SAR. Besides, let us mention space-time domain methods, such 
as the 3D Backprojection (BP) algorithm. Also known as Time-Domain Correlation, this tomographic 
reconstruction is achieved by coherent summation and stems from the projection-slice theorem; see for 
example Q. 

Back to the far-field condition of PFA, it can be relaxed even if it must be stressed that 3D near-field 
radar imaging remains computationally expensive Q, Q. Let us mention for example Q that develops 
a near-field focusing function that accounts for wavefront curvature and propagation loss; it is employed 
in Q for RCS analysis of large trees inside an anechoic chamber. Another condition, generally implicitly 
required in radar imaging, is the so-called ’’small bandwidth small angle” 0 This assumption means that 
the data is collected for a bounded excursion of the observation angles and frequency where the nature 
of the wave-object interaction is unchanged. It is closely related to the validity domain of the standard 
’’multiple scatterer model” where the object is represented by a collection of coherently illuminated point 
scatterers 0. the properties of which do not vary in the limited frequency and angle excursion band. 

In this article, we present a fast and efficient 3D non conventional radar imaging technique for identifica¬ 
tion and characterization of radar reflectivity/backscattering components of complex objects, in the context 
of RCS or signature analysis. It goes far beyond the previously proposed ad hoc localization processing 
technique, the limitations of which are described in Q. Unlike Q, our spherical measurement set-up, 
designed and developed at CEA, is especially appropriate for measuring small targets of lower RCS. The 
specific feature is that the collected scattered field data is made of polarization-diverse measurements: the 
electric field varies both at emission and reception during the acquisition, forming concentric circles with 
a singularity in the main direction. In a way, the issue is related to polarimetric imaging and to what is 
often called imaging radar polarimeter p0| |. SAR polarimetry is a widely used technique for measuring 
and identifying polarimetric properties of scatterers [ fTT| , [ fT^ . Many works have been achieved in HR 
polarimetric characterization or decomposition, especially for automatic target recognition p^ , in, in 
order to determine the full scatterer matrix of a scatterer and get a relevant information on the shape, 
orientation or dielectric properties. Here, our singular set-up implies that all the polarimetric information 
is mixed. To overcome the issue, an original approach is proposed and developed. It is based on a suitable 
extension of the multiple scattering model that leads to a high dimensional linear direct problem. This 
inverse problem is solved by a fast regularization algorithm that manages to determine at a time three 
huge 3D scatterer maps. 

The article is organized as follows. First, in section |II| a general description of the polarization-diverse 
measurement problem is made. In section |I^ is dedicated to the 3D radar imaging approach. After a 
brief introduction to classical radar imaging and the multiple scatterer model, an extension is developed. 
Inspired by HR polarimetric characterization, it leads to the expression of the direct model and to the 
proposition of a fast regularized inversion method. In section |IV| the 3D spherical RCS set-up is presented: 
it is able to accurately acquire polarization-diverse measurements. After a few details on the associated 
specific 3D radar image processing, results are presented and discussed for various mock-ups. Finally, 
conclusions are summarized in section |Vl 
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II. Problem formulation 
Standard RCS Synthetic Aperture Radar 

A monostatic radar illuminates an object with a quasi-planar monochromatic continuous wave (CW) of 
given frequency /. The incident electric field is vector of complex amplitude E^. The object backscatters 
a CW with the same frequency. The scattered electric field is vector E'^, of complex amplitude E'^. The 
complex scattering coefficient S quantifies the whole object-EM wave interaction; it indicates the wave 
change in amplitude and phase. Note that it is directly linked to the Radar Cross Section. They can be 
defined in far field condition by (see p3] | for more details): 

E‘S' 

S = lim 2^R—, RCS = |<SP (1) 

R —^cx) ej 

where R is the radar distance at which scattering is observed. The complex scattering coefficient S can be 
measured with an appropriate instrumentation system (antenna, network analyzers, etc.) and a calibration 
process. The measurements can be repeated for different wave frequencies, which is usually called 
’’Stepped Frequency Continuous Wave” (SFCW) acquisition mode. They can be repeated also for different 
incidence angles, by rotating the object or/and the antenna. Finally, it leads to a sequence of M measured 
complex scattering coefficients - - Considering all these successive measurements, it is 

generally assumed in HR radar imaging that the wave-object interaction nature is unaltered. In particular, 
the polarization configuration remains stable, i.e. the electric field orientation relatively to the target is 
unchanged. 



Fig. 1. Polarization-diverse radar imaging acquisition at successive wave frequencies {/i, /2, •' • ? /m} 


Polarization-diverse measurements 

Far from above standard RCS SAR, we consider afterwards the singular but relevant problem of a 
monostatic acquisition composed of M polarization-diverse measurements. It is represented in figure [T| 
Our goal is to analysis the scattering from Oz viewpoint. In this singular set-up, each measurement is 
associated with a different linear polarization, i.e. a different electric field direction. Let us focus on 
the i-th measurement at wave frequency /j. It is described by the wave vector kj, i.e. its direction (unit 
vector kj) and its modulus or frequency wavenumber k* = 27r/j/c (where c denotes the light speed), by the 
polarization direction at emission (unit Jones vector ef) and by the polarization direction at reception (unit 
Jones vector ef). The incident electric field E[ is collinear to ef while kj A ef indicates the direction of 
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the incident magnetic field. Finally, the polarization-diverse measurement is defined, i.e. both microwave 
emission and reception, by the successive frequencies {/i, / 2 , • • • > /m} and the successive corresponding 
directional triplets {(ki, ef, ef), (k 2 , ef, ef), • • • , (kM, e^, e^)}. 

Consequently, the polarization-diverse measurement 3D radar imaging problem consists in determining 
3D scatterer maps from Oz radar viewpoint, based on this sequence of polarization-dependent measure¬ 
ments. Compared to classical radar imaging, the main issue is linked with the variation of polarization, i.e. 
the rotation of both the electric and magnetic fields. Indeed, the classical multiple isotropic scatterer model 
on which radar imaging usually relies is no longer valid. Whereas the backscatter data is usually recorded 
separately for each polarization at emission and reception (i.e. HH, VV and HV), here all the information 
is mixed up. Besides, the huge number of unknown quantities needs to be processed efficiently. 


III. Fast 3D radar imaging with polarization-diverse measurements 


A. Overview of classical 3D HR radar imaging and the multiple scatterer model 

Let us briefly introduce classical HR radar imaging Q and first the so-called ubiquitous ’’multiple 
scatterer” model on which it is based. The multiple scatterer (MS) model lies on the interpretation of 
the target-electromagnetic wave interaction by a collection of coherent illuminated and localized point 
scatterers Related to the high frequency behavior of EM scattering from geometrically complex bodies, 
the MS model can be derived from Maxwell’s equations by high frequency approximations Q. It must 
be stressed that it requires several restrictive conditions and does not explicitly account for various EM 
phenomena: multiple scattering, creeping waves, etc. Even so, the MS model is practically and intensively 
used in radar interpretation and analysis, far beyond the validity of its assumptions. It works pretty well 
in the far field as long as there is a critical look at its output where possible artifacts may occur 0, 
mainly related to the above mentioned EM phenomena. 

MS model: In a monostatic radar context with far field/planar monochromatic wave (with wave vector 
kj), the MS model considers a target made of N elementary isotropic scatterers. They are located at points 
(n = 1, • • • ,N), of coordinates {xn,yn,Zn) in the target coordinate system (Oicyz). Each scatterer is 
defined by its complex scattering coefficient or power (n = 1, • • • , N), that quantifies the relationship 
between the incident and scattered field amplitudes (E^ and resulting from the elementary wave- 
scatterer interaction. Considering a large enough radar-target distance R, it is possible to derive from 
0 : 


E: 


1 


-{Sr, 


3 


') -E^ 


( 2 ) 


2^/7:\R ■ Ur + Tni 

where Ur is the unit vector corresponding to radar-(9 direction in ((9xyz), so that • U/j -f r„| is the 
distance between the radar and the current scatterer n. In the phase shift, the term 2ki • accounts for 
the back and forth delays. 

In the MS model, E'^ = J2n=i received echo or scattered signal is ideally formulated as a 

coherent superposition of elementary echoes 0, neglecting multiple scattering terms. That leads to the 
following scattering coefficient: 


N 


s,. 


T?S 

= 2v^«|7 = E 

n=l 


R 




Since the ratio 


R 


IR-UH+r^l luH+rn/RI 


R-UR + Tnl 
the MS model reduces to: 


(3) 


N 


<5.=E' 


3 


(4) 


n=l 


The N complex scattering coefficients associated to the multiple scatterers determine the system 
response. They do not depend on the wave vector kj, i.e. its amplitude and direction. 


'it is assumed that the incident scattered field amplitude does not depend on n. 






5 


Typical application: Consider the typical acquisition of Fig. that can be encountered in spherical 
3D radar imaging. It requires a sequence of M viewpoints around Oz. Considering that the i-th wave 
direction is defined by the standard spherical angles (r?*, $ 1 ), the successive wave vector directions kj are 
given by: 

kj = [— cos Dj sin $j — sin t^j — cos t?* cos $*]* 


It must be noticed that in such an acquisition the electric field varies slightly as long as the variations 
of (t9j, $j) are limited. This is also called the ’’small angle” condition. Moreover, the electric field stays in 
the same plane, as shown in Fig. and the polarization remains more or less constant. Another condition 
is required: it is called the ’’small bandwidth” condition. It implies that the target-wave interaction does 
not significantly vary within the limited frequency variation domain. If these conditions are fulfilled, each 
complex scattering coefficient is considered constant for all the measurements. The MS model Q can 
be formulated for each measurement i (for i = 1, - • • , M) as: 


N 


■Si = E' 


n=l 


^ ^2jki-(smi^i cos^iXn+sini^i sin^iyn+cosi^iZn) 



Fig. 2. Classical spherical radar imaging with ef = ef = 

Remark that the MS model can be formulated as: 

+ e (5) 

where X(r) = — r^) is the spatial density of scatterers and (i(-) is the Dirac function. The 

additive term e accounts for noise and uncertainty, due to interfering echoes generated by the environment 
and MS model limitations. It is generally considered as white and Gaussian (centered, complex and 
circular). 

Radar imaging: Radar imaging processing is an inverse problem. For a given radar viewpoint (e.g. 
Oz), it consists of the determination the complex 3D image or map X, i.e. the complex amplitudes Sj at 
points Tj. The above mentioned ’’small bandwidth small angle” condition is assumed. It must be stressed 
that it requires that the wave polarization remains constant so that the target-wave interaction does not 
vary from one measurement to another. 

The Polar Format Algorithm (PFA), formerly introduced, is based on Q which is the discrete Fourier 
transform (DFT) of the sought image X. The computation of X is achieved via the inverse Fourier transform 
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of the measured hologram, Sj (for j = 1, - , M). In actual practice, the limited and discrete acquisition in 

angles and in frequencies must be taken into account. For a regular acquisition in angle and frequency, an 
interpolating pre-processing step of regridding is performed before an inverse 3D Fast Fourier transform. 


B. Extension to the multiple scatterer model 

Considering the polarization-diverse acquisition of Fig. the MS model is no longer valid and as a 
consequence conventional radar imaging processing can not be applied. Actually, the complex scattering 
coefficients (n = 1, • • • .N) can no longer be considered constant. The electromagnetic interaction 
strongly depends on the electric field direction; consider for example the scattering with the edges of an 
object. To overcome this issue, we next develop an adaptation of the MS model in this varying polarization 
context that will latter support the imaging inversion process. 

General formulation: Remind that the purpose is to achieve radar imaging analysis from Oz radar 
viewpoint. Let us consider the i-th acquisition of Fig. [T] around Oz. Then, the basic idea is to consider 
that each elementary isotropic scatterer is associated with not only a complex scattering coefficient but 
with a full polarization scattering matrix S^: 


S„ = 


XX 

n 

,xy 

n 


xy 

n 

yy 

n 


( 6 ) 


More exactly known as the radar cross section matrix p^ , S„ is defined in the target reference frame 
Oxy, associated with an illumination direction collinear to vector z. It quantifies the amplitude and 
polarization of the scattered wave for an arbitrary polarization of the incident wave. Note that 
due to the reciprocity theorem for a monostatic case (T^, so that the scattering matrix S„ is completely 
defined by the three complex coefficients: s^, and s^. Next, we derive an extension to the MS model. 
For each illumination, it introduces the appropriate combination of scattering matrix terms. 

Introducing polarization, the former general definition ([T]) of the complex scattering coefficient can be 
reformulated by: 

rj^S 

S = lim 2\/t:R— (7) 


R —^(X) 


'E^ 


where Gr is a unit vector aligned along the electric polarization at reception and the scalar product 
= Gr- E'^ corresponds to the transmitted electric field. Refer to for more details. 

The polarization configuration of the n-th scatterer is represented in Fig. For the z-th acquisition, the 
scattering field vector resulting from the interaction between the incident wave and the n-th elementary 
isotropic scatterer is given in (Ox'y') by (for a large target-radar distance R): 


E 


s 

n(x',yO 


g-2jki-r„ 

2^\R • Uii -F r„ 

g-2jki-rn 

2^\R • u/i + r„ 


(SU*) • E: 


'i (x',y')J 



where S^(z) is the scattering matrix in the reference frame Ox'y', normal to the current direction of 
propagation kj. Note that there is no component of E^ and Ef in the direction of kj. 

Let denote ef the unit Jones vector at emission and ef the unit Jones vector associated with the 
receiving linear polarized antenna. The radar received voltage is proportional to the transmitted electric 
field: Tf = [efj., so that: 




rpo _ 

^ n 


2y/w\R • Ui^ + 

" 2^/^T\R ■ Ur r„ 


ef ■ x' 

_ ef • y' _ 

' • s;(^) ■ 

1 1 

H H 

1_1 

■ ef • x' ■ 
ef • y' 

* • ■ 

■ ef • x' ■ 
_ ef • y' _ 


It comes to: 
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Fig. 3. Interaction of an EM wave with a point scatterer 


Similarly to the standard MS model, the received echo can be simply formulated as a coherent 
superposition of the elementary echoes: T;^. Thus, the complex scattering coefficient scattering 

coefficient is: 


rpS rri, 

S, = 2V^R-^ = 


( 8 ) 


In the far field {Vn/R ~ 0), it is then straightforward to derive the following extended MS model for 
polarization-diverse measurement, incorporating noise: 


N 




i*(i) • e 


+ 


n=l 


with slii) = [ef • x' ef • y'] • S;(i) 


ef • x' 
ef -y' 


(9) 


( 10 ) 


As previously with the MS model about the scattering coefficients s^, a ’’small bandwidth small angle” 
will be further assumed about the scattering matrices. Due to a limited acquisition domain, in angle 
(around Oz) and frequency, the scattering matrices are supposed to be nearly constant: S„ for 

i = 1, • • • , M. Former expression ([T^ leads to: 


Sn(0 ~ [ef • x' ef • y'] • S 


ef • x' 

ef-y' 


( 11 ) 


Basically, this expression provides the complex scattering coefficients or powers s* of the MS model that 
are adapted to the i-th configuration. Note again that it depends on the relative target attitude towards 
the emitting and receiving antenna. In classical radar imaging, the directions of ef and ef remain stable, 
selecting constantly a term of the scattering matrix or a fixed linear combination of them. In this 
varying polarization context, they are all ’’mixed”: and ( fTT| ) show that each acquisition is to select a 

different combination of the scattering matrix S„. Next, the extended MS model is detailed for a specific 
set-up. 

Specific form for concentric polarization-diverse measurements: Let us consider the concentric acqui¬ 
sition, with the successive wave and electric field vectors represented in Fig. It stems from a spherical 
3D RCS set-up that will be further detailed in section |IV| Firstly, note the following sort of singularity: 
the vectors spin around axis z and there are various acquisitions from wave vectors colinear to z while 
electric fields differ. Note that it induces different scattered fields, it goes against the classical MS model 
with scalars s„. 
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Fig. 4. Concentric Polarization-Diverse Radar Imaging Acquisition 


The i-th measurement is given by for i = I,-- - ,M. The azimuth 6i and the roll cpi 

correspond to two rotation angles (the roll rotation is defined around z, see section |IV] for details), that 
defines the directions of the wave vector direction kj as well as the emitting/receiving linear polarization: 



— sin 6i cos (fi 

— sin 9i sin ipi 
— cos 9i 


(x,y,z) 


— cos 9i cos Pi 


— sin Lpi 

— cos 9i sin Pi 

or 

cos ifi 

sin 9i 

(x,y,z) 

0 




depending on the emitting and receiving linear polarized mode of the antenna (respectively H or V), and 
so for ef. 

Afterwards, we determine the specific forms of (ef • x', ef • y') and (ef • x', ef • y') in order to express 
s*(i) for concentric acquisition. First of all, let us define the intermediate coordinate system (xQ,yQ,ZQ): 


yo 


-k* 

z'o A z 

II ^0 A z II 

Yo A z'o = 


— sin (pi 
cos ipi 

0 

cos 9i cos ipi 
cos 9i sin ipi 
— sin 9i 


(x,y,z) 


(x,y,z) 


( 12 ) 

(13) 

(14) 


The unit Jones vector at emission ef, and so at reception ef, can be expressed in (xg^yo) = 

(V mode). The antenna coordinate system (x',y',z') is defined 

by: z' = Zg = —kj, x' is collinear to the projection of x on the orthogonal plane (xg,yg) and y' = z' Ax'. 


cos 9i cos (/PjXg — sin (piYQ 
\J cos^ 9i cos^ Lpi + sin^ (/?* 
sin 9 ?jXo + cos 9i cos 
\J cos^ 9i cos^ Pi + sin^ pi 


(15) 


( 16 ) 
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For the H mode, it results in the following expressions for ef ■ x' and ef ■ y' (and so ef): 




e^-x' = 


e^-y' 


— cos 9i cos ipi 


\l cos^ di cos^ Lpi + sin^ Lpi 

— sin Lpi 

a/cos^ 9i cos^ ifi + sin^ ipi 


(17) 

(18) 


Let us sum up the complex scattering coefficient in this concentric polarization diverse acquisition. For 
the i-th polarization acquisition, the complex coefficient s*(i) is given by: 

1) HH acquisition mode 


4(0 


ef-x' 


ef • x' 

4"^ • y' 


^ [cos^ 9i cos^ ipi ■ + sin^ ipi ■ s 

+ cos 9i sin 2(pi ■ 


yy 

n 


(19) 


2) VV acquisition mode 


3) HV acquisition mode 


~ ^[sin^ ipi • + cos^ 9i cos^ ipi • s)7 

— cos 9i sin 2ipi ■ 


*/■^ 1 r .sin 29 ? sin 29 ? 

SnW ~ -j^[cos9^—■ -cos9^—■ 

— (cos^ 9 cos^ ip — sin^ p) ■ 


( 20 ) 


( 21 ) 


with: JCi = cos^ 9i cos^ pi + sin^ pi. 

Afterwards, we consider the whole measurements, i.e. the M observed complex scattering coefficients, 
in the associated acquisition conditions defined by the successive frequencies fi, azimuth 9i and roll pi 
for z G 1, • • • , M. They are stacked in the following complex vector S: 


S=[Si S2 ••• <SMf 


( 22 ) 


C. Derived direct matrix model 

We now derive the observation model (forward model) from the previous extended MS model. Let 
us stress again that the ’’small bandwidth small angle” is assumed: the scattering matrices are supposed 
to be nearly constant for each elementary scatterer. The unknown is made up of the scattering matrices 
associated to the N discretized scatter points r„ (n = 1, • • • , A) of a 3D grid, the coordinates of which 
are {xn-,yn-,Zn) in the target coordinate system (Oxyz). The scattering matrices can be separated in 
three 3D maps Sxx, Syy et Sxy (in the target reference frame), respectively for HH, VV and HV antenna 
polarization mode (defined at the initial reference rotation, for 9 = and p = 0°). In vector form, 
they can be expressed respectively W: Sxx = [-s^ '^ 2 ’^ ' ‘ ‘ ^yy = Wi ^ 2 ^ ''' 

Sxy = [si^ • • • s^]*. From (%, the complex scattering coefficients Si are clearly related to the 

Fourier transform of the s* coefficients. Furthermore, it is obvious from ( [T9] ) and ( |20| ) that the scattering 
coefficients are linear combinations of components Sk (k G {xx, yy, xy}). Therefore, the direct model (j^ 
can be rewritten as (with noise vector e): 

S = 'y ^ TWkF^k + ^ 

k£{xx,yy,xy} 


( 23 ) 
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where F is a DFT (Discrete Fourier Transform), that can be computed by FFT (Fast Fourier Transform) al¬ 
gorithm. Finally, T is the truncation matrix related to the position of the observed points in Fourier domain. 
The three matrices Wk are diagonal weight matrices that affect Fourier coefficients (k G (xx, yy, xy}): 

Wfc = diag(wfc(l),Wfc(2), • • • ,Wk{M)) (24) 

Depending on the polarization acquisition mode (HH, VV or HV), the weights Wk are given respectively 
by ( fT9| ), ( p0| ) and 0. They define the elements or combination of elements of the N MS scattering 
matrices (i.e. s^, or , for n = 1, • • • ,N) which determine the scattering coefficient observation <Sj. 
When the acquisition mode is HH, the weights are given by (for z = 1, • • • , M): 

{ Wxx(«) = COS^ 6i COS^ 

Wyy(z) = sin^ (25) 

Wxy(^) = cos^j sin2(/7j//Cj 

They are shown in Fig. for various acquisition angles 6 and (p. Typical acquisition angles {6 < 30° 
and 0° < (/? < 360°) are represented by black points. Notice that the weights Wxx, Wyy and Wxy depend 
mainly on the value of p. For p is close to 0°, only the terms matter. Conversely, when p close to 
±90°, only the matter and for that reason are observable. Notice that it is true as long as 9 < 40°. For 
larger 9, the scattering terms no longer affect the coherent superposition echo. Anyway, the ’’small 
bandwidth small angle” assumption of constant will be violated for common complex targets. 



theta (”) 



theta (°) 



20 40 60 

theta (°) 


Fig. 5. Weigths Wxx, Wyy and |wxy| in function of 0 and if, with typical acquisition angles (black points). 


Finally, merging the 3 maps in an unique column vector s = [sxx; Syy; Sxy]> the final direct model reads: 

S — As ± £ (26) 


where A is a large block matrix of size M x ZN: 


A = [ TWxxF I TWyyF I TWxyF ] 


(27) 


The observation matrix A defines the deterministic part of the relation between the observations and the 
unknown state vector made of the three vectorized 3D map. It must be noticed that its dimensions can 
be huge, e.g., N k, 500000 and M = 15 • 10®. And yet, the model (^) is still linear (linear transforms 
and additive noise). Moreover, it relies on simple and fast transform (FFT, weight, truncation). Again, 
remark that a close parallel can be made with polarimetric SAR p0| |, p^ , with a linear and additive 
noise model and a polarimetric measurement matrix which, similarly to A, is determined by the successive 
polarizations at emission and reception. 



















































11 


D. Regularized inversion and Minimum Norm Least Squares 

Regarding the inverse problem, a co m mon approach relies on a discrepancy between measurements 
and model outputs. A standard discrepancy is based on a quadratic norm and yields the so called Least 
Squares (LS) criterion: 

= ||5-Asf . (28) 

A potential solution minimizes Ji^s, i-e., it is the one that best fits the data: 

Sls = arg min Jls(«) • (29) 

seC^ 

The minimizer Sls can be found by setting the gradient of to zero and it leads to a linear system: 

A^As = A^S. (30) 

Unfortunately, it cannot be solved since A^ A is rank-deficient (the number of unknowns is larger than 
the number of measurements) and, as a consequence, an infinity of objects minimizes (and nullifies) the 
LS criterion. In other words, an infinity of backscattered maps is exactly consistent with the data. Among 
these solutions, a possible simple approach selects the one with minimum norm: 

|2 

(31) 


5mnt,s - 


arg mm | 
se€^ 

s.t. «S — As = 0 


i.e., the Minimum Norm Least Squares (MNLS) solution. There is a very large literature on various 
methods to deal with this kind of situations and build other solutions. The interested reader can, for 
instance, refer to books on inverse problems [fT7||-[|2T]|. Coming back to our simple MNLS solution, it 
can be shown that 

= At(AAt)-i5. (32) 

A possible proof is based on Lagrange multipliers and it is given in Appendix. The estimated map Smnls 
is a linear transform of the data, nevertheless the matrix AA^ is huge and its inverse cannot be computed. 
Fortunately, it can be inverted analytically: given the specificity of ( |27l ), we have 

AA^ = T {Wl + Wfy + Wj) T* 


that is a diagonal matrix and hence easily invertible. Finally, relation ( [32| ) becomes 

^ FtWxxT* ^ 


^MNTA - 


FtWyyT* 

FtWxyT* 


[T{W, 


+ w‘ +W‘ 

XX ' yy xy. 


T*] 


that jointly determines the three 3D backscattered maps associated with the xx, yy and xy polarizations. 
In addition, each individual map obtained separately by: 

= F'WtT [T (Wj + + wD r‘] 5 

from the data set S. It can be seen that they are obtained by IFFT of the zero-padded weighted measure¬ 
ments. The involved weights are: 

Wk{m) 


=- ! ^2 , - T^TL - TNd. 

+ Wyy[m)^ + w^y[my 

based on direct weight ( |25| ), that depends on the polarization mode and the acquisition angles (azimuth 
and roll). 

The resulting fast 3D radar imaging algorithm is composed of two steps as described in Fig. The first 
step consists in processing the data by applying the weights of Eq. ( |^ , for each polarization xx, yy and xy, 
and for each associated weights. Afterwards, in the second step, three 3D backscattered maps, associated 
with XX, yy and xy (in the target reference frame), are separately reconstructed by a 3D PFA method 
[ [22| |. It includes regridding and inverse 3D IFFT steps. It is finally a very efficiently implementation, 
based on weight, zero-padding and IFFT. 
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Measurement angles 
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Fig. 6. Fast regularized inversion 


IV. Application 

Indoor 3D spherical near field facilities are dedicated to various tasks. Beyond scattering analysis with 
2D or 3D ISAR imaging, they are encountered for RCS measurement, microwave absorber measurement, 
nondestructive testing, material characterization, antenna measurement and diagnosis, research in near¬ 
field measurement techniques, etc. These facilities are generally based on a mobile arch that provides 
various illumination viewpoints. Most of them are dedicated to antenna measurements; let us mention for 
example the commercial systems ”G-DualScan” and ”7008-90”, respectively developed by Satimo and 
NSI. Others are specialized in RCS measurement. See, for instance, the system of the French research 
institute Fresnel at the Microwave Common Resources Centre (MCRC), equipped with five positioners 


that allows spherical diffraction patterns p3] |, as well as the bistatic anechoic chamber BIANCHA p4 | 
of the Italian institute INTA. It provides a bistatic, spherical field scanner based on a dual-axis azimuth 
turntable and two elevated scanning arms. Related to these 3D spherical facilities, our 3D spherical RCS 
set-up, designed and developed at CEA, is specifically dedicated to RCS analysis of small targets of lower 
RCS. After a short description of the facility and its microwave instrumentation, we present 3D radar 
imaging results on various mock-ups. 


A. The 3D RCS CEA facility 

The concentric polarization-diverse measurements of previous section |III-B| are achieved by the exper¬ 
imental layout of Fig. Dedicated to the RCS characterization of small targets (typically < 1 ft), it is 
composed of a 4 meters radius motorized rotating arch (around horizontal axis) holding the measurement 
antennas while the target is located on a polystyrene mast mounted on a rotating positioning system 
(around vertical axis). The combination of the two rotation capabilities leads to a full 3D near field 
monostatic RCS characterization. 

The microwave instrumentation is made of two bipolarization monostatic RF transmitting and receiving 
antennas that are driven by a fast network analyzer. A phased array antenna is used for the 0.8 -1.8 
GHz frequency band. It is optimized in order to reduce spurious signals originated from interactions with 
the arch metallic structure. For higher frequencies from 2 to 12 GHz, a wideband standard gain horn, 
equipped with a lens, is preferred (see Fig. |^. The target under test is located vertically on a polystyrene 
mast below the measurement antenna. The target rotation around vertical axis is determined by the roll 
angle (/?. The arch rotation around horizontal axis is determined by the azimuth angle 9. A full set of 
data S is therefore composed of successive frequency, azimuth, and roll sweeps. The arch is built of 
aluminum, weights about 200 kg, measures approximately 5 m x 3 m, and is designed to minimize 
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Fig. 7. Near-field 3D RCS measurement facility 


potential mechanical deformations and to hold the measurement antenna over a total travel range of 
±100°. The 6 rotation is achieved using a direct drive positioner including a brushless motor with a high 
accuracy encoder directly coupled to the motor for accuracy and repeatability concerns. The ip rotation 
uses another positioner also including a brushless motor with a high accuracy encoder directly coupled 
to the motor. The polystyrene mast p5] |, [ [26| supporting the target under test is located on this rotating 
positioning system (vertical axis). Let us emphasize the precision of the whole positioning system. The 
two major errors are the residual radial error and the repositioning (repeatability) error of the arch. The 
first one is about ±0.001° and is equivalent to a maximum error of ±0.75 mm on the vertical axis. The 
second one is about ±0.15 mm maximum, once again on the vertical axis. Refer to p7] | for details on 
the laser tracker characterization. 

Concerning the measurement process, it is composed of successive stages. First, a calibration by 
substitution is performed: it consists in replacing the test target with a calibration standard whose 
echo is well known in order to determine the inverse transfer function /C of the entire RCS measurement 
system. This enables to convert indirectly the received transmitted electric field Tout to a complex scattering 
coefficient value: S = fC.Tout- Basically, the substitution method establishes a phase reference relatively to 
the rotation center and normalizes out the dispersive frequency response of the RCS measurement system. 
Note that for large objects, it can be enhanced by a multi-calibration approach p8| | which overcomes most 
of near-field effects, such as the power decrease or the EM wave sphericity. The second stage corresponds 
to background subtraction [ fT^ . It is commonly introduced when the clutter is too high, compared to the 
test target and calibration signals. It removes coherently the background echoes. Eventually, an adaptive 
range filter is applied to eliminate the residual interferences, e.g. interactions with walls and floor that 
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can affect the useful signal. 

B. Results 

Let us first illustrate our 3D radar imaging method by an application to the following synthetic RCS 
data. We consider the metallic ogival-shaped object from the EM benchmark [ |29| |, with polarization- 
diverse measurements at 9 e [0° : 2° : 20°], (/p G [0° : 5° : 360°[ and / G [1 GHz : 10 MHz : 3 GHz]. The 
synthetic data is computed by an efficient parallelized harmonic Maxwell solveiQ Fig. shows the 3D 
scatterer map corresponding to xx polarization. The 3D map corresponds to the indicated above 

radar viewpoint, along Oz. Here, the 3D volume is represented by several transparent isosurfaces, that 
are linked to different scatterer levels. In this basic situation, the main EM scatterers are located: the 
diffraction on the above tip and, at a higher level, the diffraction on the below tip. It comes with sidelobes 
that derives from the limited acquisition sampling. 






-41 dB 
-51 dB 
-61 dB 


Fig. 8. Ogival shape 3D radar map ^^nls = 64 x 64 x 1024, M ps 2 • 10®] 

Next, it is shown how the 3D radar imaging approach can support RCS analysis in more complex 
situations from real polarization-diverse measurement data. Various targets have been characterized in 
order to evaluate the 3D radar imaging method: a metallic cone with patches, a metallic arrow and a 
glider. Each of them enables to test the capacity to deal with various target-wave interaction phenomena: 
localized scatterers with the metallic cone, multiple diffractions as well as different polarization responses 
with the metallic arrow and near field effects due to the large wingspread of the glider mock-up (regarding 
to the wave number). All the scattering measurements have been acquired for a 360° roll sweep (cp) and 
a ±20° azimuth sweep (9), from 8.2 to 12.4 GHz frequency. 

1) Localized scatterers: Fig. shows the 3D radar map polarization xx of a metallic cone 

(height: 60cm) where 3 metallic patches (lcmx4cm) have been glued to points z = 250 mm (—45° roll), 
z = 400 mm (135° roll) and z = 112 mm (270° roll). For convenience, the 3D volume is represented by 
several isosurfaces, superimposed upon a volume section relatively to the target shape. The main scatterers 
are perfectly located: the tip, the diffraction of the rear edge and each metallic patch. Notice that, since 
the frequency band is here larger than for the previous ogival-shaped object, the resolution is improved 
along Oz. The scatterers corresponding to the 3 metallic patches are correctly located. 


^It combines a volume finite element method and integral equation technique, taking benefit from the axisymmetrical geometry of the 
shape (30) 
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Fig. 9. Metallic cone with 3 metallic patches (left) - 3D radar map = 256 x 256 x 512,10^] (right) 


It is confirmed in the above view of Fig. [T^ where the image can be compared to the photo, as well as 
in the exploded image of Fig. There, the true patch locations are represented by black pellets. Notice 
that although the creeping waves are not shown in Fig. they would appear, under the cone, in other 
representations. 



Fig. 10. Exploded image of the 3D radar map Smnls [N = 256 x 256 x 512, M w 10®] 
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Fig. 11. Above view of the metallic cone with 3 metallic patches 



Next, the metallic patches are replaced by 3 RAM (Radar Absorbing Material) rectangular patches 
(lcmx4cm) which have been glued to the same locations. Results are presented in Fig.[T^ If the scatterers 
associated to the 3 RAM patches are correctly located, they are now far more diffuse; it is then more 
difficult to evaluate the location accuracy of the 3D imaging method. This phenomenon can be explained 
by the fact that RAM has a spread frequency response, which does not allow to obtain a very localized 
scatterer. 

2) Multiple diffraction & polarization dependence: The second target is a metallic arrow of 472 mm 
long and 150 mm large (it belongs to Airbus Group Innovation, see [ |3T| | for details). It is represented 
in Fig. [T^ With such a target, diffraction edges are known to differ according to the polarization of the 
incident wave. Moreover, multiple diffractions are likely to appear between the back of the arrow and its 
base. 

Fig. [T^ provides a comparison of 3D radar maps of the arrow in polarization xx. Radar imaging from 
real data is compared to radar imaging from simulated data, computed with a parallelized 3D EM solver 
(see for details). It can be checked that both maps are very close one another. The sensitivity to wave 
polarization is emphasized in Fig. [T^ and [T^ The xx polarization map of Fig. is again represented in 
the shape section of Fig. [T^ Depending on the polarization (i.e. xx, yy or xy), the diffraction on the edges 
is more or less important, as well as the multiple reflections of the electromagnetic wave between the rear 
of the arrow and its base. Located below the target, they are much more important in and s^nls. due 
to the higher resonance with the base. Let us stress that the 3D Radar imaging enables to localize not 
only 3D scatterers associated to diffraction and reflexion but also scatterers that depolarize the incident 
wave, like the corners of the target. So the simultaneous determination of the three 3D scatterer maps is 
particularly helpful for RCS analysis: it permits to localize and characterize each 3D backscatterer. 
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Fig. 13. Metallic cone with 3 RAM patches - 3D radar map Smnls (shape section) [Ai = 256 x 256 x 512, M « 10®] 



y (m) X (m) y (m) 


x(m) 


Fig. 14. 3D radar image (Smnls) of tho arrow, from real data (left) and simulated data (right) [A^ = 40 x 40 x 1146, M k. A - 10®] 


3) Near field effects: They can be shown with the full aluminum glider-shaped mock-up of Fig. |7| and 
its large wingspread. Its main dimensions are approximately 3 ft x 3 ft x 8 in (see for details). The 
glider mock-up is located vertically on the polystyrene mast. 

In Fig. [T^ and 3D scatterer maps are shown depending on polarization, respectively xx and yy. 
Even if the maps are a bit more complex with several multiple scattering artifacts, the main scatterers of 
the target are apparent: the nose, the fuselage and the tail. However, it can be clearly seen that it does 
not exhibit the expected high specular behavior of the trailing edge of the front wing glider. Contrary to 
far field, the backscatterrer is not located uniformly along the entire length of the wing. It mainly results 
from aliasing, due to angular undersampling both in 0 and cp, in the context of the extended wingspread in 
the transversal dimension. It also partly results from the spherical RCS data and the anisotropic radiation 
pattern of the antenna These defaults can be corrected by improving the angular sampling and by 
introducing a 3D near field/far field correction in the radar imaging algorithm. 


V. Conclusion 

A 3D radar imaging technique has been presented. It is able to process a collected scattered field data 
made of polarization-diverse measurements, where the electric field direction varies during the backscatter 
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Fig. 15. 3D radar image of the arrow (from real data): sJJnls (left) and (right) [At = 40 x 40 x 1146, M k. A - 10®] 



Fig. 16. 3D radar image of the arrow: Smnls (shape section) [At = 40 x 40 x 1146, M « 4 • 10®] 


data acquisition. It is based on fast and efficient regularized inversion that reconstruct three huge 3-D 
scatterer maps at a time. The approach rests on a simple but original extension of the standard multiple 
seatterer point model, elosely related to HR polarimetrie eharaeterization. It is applied sueeessfully to 
synthetic and real data, colleeted from a 3D spherical experimental layout dedieated to aeeurate RCS 
characterization. 

To go further, the resolution could be enhanced by introducing non-quadratic approaches (L2-L1) p4| , 
p5] | or eonsidering none differentiable criteria. From p6| , eonstraints could be introdueed, e.g. from the 
known shape, in order to improve the seatter inference and remove ambiguities Q. Besides, similarly to 
it eould be possible to take into account the spherieal near field illumination of the target as 
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well as the radiation pattern of the antenna. 
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Appendix 

Consider the general quadratic problem in with P linear equality constraints: 


X = 


argmin ||a?||Q 

a;eCJV 

s.t. c — Ax = 0 


(34) 


where 
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• Q is a X positive-definite matrix that defines the criterion to be minimized 

• while the constraint is decribed through the vector c G and the P x N matrix A. This matrix is 
assumed to be of full rank. 

A standard solution to such a problem relies on Lagrange theory: Lagrange multipliers, duality and saddle 
point, as follows. The Lagrangian of the problem ( |^ reads: 

C{x,u) = x^Qx + u\c—Ax) (35) 

where w e is the Lagrange multiplier, also referred to as the dual variable {x G is referred to as 
the primal variable). 

Let minimize C with respect to x (for a fixed u) by setting to zero the gradient of the Lagrangian: 

dC 

(x.u) = 2Qx — A^u = 0 
ox 

The solution is directly given by: 

X = ^Q~^A^u. (36) 

Note that the Hessian is equal to 2Q and hence is strictly positive. Consequently, x corresponds to a 
minimum. Then, the introduction of its value in the expression of £ given by ( [35] ) provides the dual 
function C: ^ 

C{u) = inf C{x, u) = £(*, u) = — v) AQ~^A}u + v)c. 

X 4 

Next, let maximize £(w) relatively to u by setting to zero also the gradient with respect to the dual 
variable: 

8C 1 

0 = 7 ^(w) = —-AQ~^A^u + c 
ou 2 

The solution is directly given by: 

u = 2{AQ-^A^)-^c. (37) 

Note that the Hessian is equal to —{AQA^)/2 and hence is strictly negative; u corresponds to a maximum. 
Let insert the expression ( [J7| ) of the dual optimizer into the expression ( |36| ) of the primal optimizer. It 
leads to: 

X = Q-^A\AQ-^A^)-^c (38) 


which is the Minimum Norm Least Squares (MNLS) solution. When Q is the identity matrix, it results 


in the expression (|32|): 


c. 


Note that, if A/^ = P, A is invertible and x = A ^c. 



